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The magnetization of iieavy lioles in III-V semiconductor quantum wells with Rashba spin-orbit 
coupling (SOC) in an external perpendicular magnetic field is theoretically studied. We concentrate 
on the effects on the magnetization induced by the system boundary, the Rashba SOC and the 
temperature. It is found that the sawtooth-like de Haas-van Alphen (dIfvA) oscillations of the 
magnetization will change dramatically in the presence of such three factors. Especially, the effects 
of the edge states and Rashba SOC on the magnetization are more evident when the magnetic 
field is more small. The oscillation center will shift when the boundary effect is considered and 
the Rashba SOC will bring beating patterns to the dHvA oscillations. These effects on the dHvA 
oscillations are preferred to be observed at low temperature. With increasing the temperature, the 
dHvA oscillations turn to be blurred and eventually disappear. 



I. INTRODUCTION 



The physics of low dimensional electron systems in 
high magnetic fields is one of the most important sub- 
jects of semiconductor physics. Recently, great atten- 
tions are taken to the magnetic properties of two dimen- 
sional (2D) systems in a strong perpendicular magnetic 
field both experimentally and theoretically P, [3, S13, 01 ■ 
Such systems are candidates for next-generation spin- 
tronic devices^, 0, [1| because of their very high electron 
mobility compared with silicon. In the absence of disor- 
der and interactions, due to the Landau energy quanti- 
zation in an external magnetic field, the magnetization is 
predicted to oscillate periodically in a sawtooth pattern 
as a function of B or filling factor v, which is the fa- 
mous de Haas- van Alphen (dHvA) oscillation (ol. fiol. [T]]| . 
The dHvA oscillation is a powerful tool used extensively 
to determine the electronic properties of bulk semicon- 
ductors and metals, such as Fermi surface geometry and 
effective masses. It has been evident now that the magne- 
tization, which turns to be zero for classical electrons, ap- 
pears due to the presence of the sample boundary , the 
spin orbital coupling (SOC) or Zeeman splitting. Once 
a perpendicular magnetic field was applied, these factors 
would change the magnetization dramatically. Firstly, 
in a finite size sample with boundaries the edge states 
play important roles in quantum transport in high mag- 
netic fields, such as quantum Hall effect or magneto- 
tunneling effect. From the semiclassical point of view 
when electrons are located near the edge of the samples, 
cyclotron motion cannot fulfill the entire cyclotron or- 
bit if the distance between the center of the cyclotron 
orbit and the edge is smaller than the radius of the cy- 
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clotron motion. In such a case, the electron orbit is 
bounced at the edge plane, and the bounced electrons 
conduct a one-dimensional motion as a whole being re- 
peatedly reflected at the edge plane, resulting a skipping 
orbit. The edge states with a skipping orbit are differ- 
ent from those in the interior of the sample where the 
cyclotron orbit is completed. Although these edge states 
are localized in the vicinity of the confining walls, the 
measurable bulk quantities are completely modified by 
the edge states, certainly including the magnetization. 
Secondly, the Rashba SOC, competed with the Landau 
splitting and the Zeeman splitting, will modulate the en- 
ergy spectrum. As a result, the Rashba SOC will have 
non-neglectable influence on the dHvA oscillations of the 
magnetization. 

However, to our knowledge there are no detailed treat- 
ments on the influence of edge states and SOC on the 
magnetization in 2D systems, except for the 2D electron 
systems much recently studied by our group p^. In this 
paper, we study systematically the thermodynamic mag- 
netization of heavy holes in III-V semiconductor quan- 
tum wells with boundaries in the presence of the Rashba 
SOC due to structure-inversion asymmetry, with an ex- 
ternal magnetic field applied perpendicularly. We will 
show that the dHvA oscillations of the magnetization 
change dramatically due to the presence of SOC and sys- 
tem boundaries. The edge states lead to the change of 
both the center and the amplitude of the sawtoothlike 
dHvA oscillations of the magnetization. The SOC mixes 
the spin-up and spin-down states of neighboring Landau 
Levels into two unequally spaced energy branches, which 
further changes the well-defined sawtoothlike dHvA os- 
cillations of the magnetization. These effects on the mag- 
netization can be observed at low temperature in exper- 
iment. With increasing the temperature, the dHvA os- 
cillations of the magnetization turn to be blurred and 
eventually disappear. 

This paper is organized as follows: Secllll focuses on 
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the quantum mechanical solution of the heavy-holes in 
III-V semiconductor quantum wells. Using numerical di- 
agonalization of the Hamiltonian in a truncated Hilbert 
space we calculate the energy spectrum. Sec lIIII gives 
the main results of this paper on the magnetization and 
makes a detailed discussion. The last section gives a short 
summary. 



II. ENERGY SPECTRUM OF THE HEAVY 
HOLES 

Now we consider a 2D hole gas system in which a 
Rashba SOC arises from the quantum well asymmetry 
in the growth (z). The Hamiltonian for heavy holes in 
III-V semiconductor quantum wells within a perpendic- 
ular magnetic field B ~ ^Biz, taking into account the 
kinetic energy and the Zeeman magnetic field, can be 
written as 



2m* 2h? ^ ^ ^ 



--gsUBBuz+V {y) , 
(1) 



cient due to structure inversion asymmetry across the 
quantum well grown along the [001]-direction chosen as 
the z-axis. For a symmetrically grown quantum well, 
the coefficient a is essentially proportional to an electric 
field ap plie d across the well and therefore experimentally 
tunablelli, [H, [H, [13, IH, d. ffs is the gyromagnetic 
factor. The last term V (y) is the lateral confining po- 
tential. For simplicity, a hard wall potential that confines 
the holes in the transverse y direction is used. 
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It is convenient to use the Laudau gauge A ~ {yB, 0, 0) 
and to write the wave function in the form 



^{x,y) = e'''^^{y) 



(3) 



with the function ip (y) expanded in the basis set of the 
infinite potential well. 



where m* being band effective mass, Il± = Ilx±illy with 
= + (e/c) Arj^ Hjj and A^j denoting the rj compo- 
nent of the momentum and vector potential, respectively. 

B=V X A. The Pauh matrices, a± = ± iay, operate Introducing the cyclotron center yp 
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on the total angular momentum states with spin pro- 
jection ±3/2 along the growth direction; in this sense 
they represent a pseudospin degree of freedom rather 
than a genuine spin 1/2 [ij. a is Rashba SOC coeffi- 



the magnetic length l(, ~ ^JhcfeB and the cyclotron 
frequency ujc = e-B/m*c, the Schrodinger equation 
{x,y) = E'ii {x,y) leads to the following equations 
for spinors: 
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sin (nt) dt, 
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3{t-t„) + 3{y~ yo) dt + 



sin (nt) dt, 



T 



sin (nt) dt, 



with e = E/hoJc, t = Try/L, to = i^yo/L, g = m*gs/2me, 
and ?; = m*alb/h'^. Taking the same method in Ref . [20| . 
we solve the above equations in a truncated Hilbert space 
disregarding the highest energy states. Typically we take 
a matrix Hamiltonian of dimension of a few hundred. We 
increase the size of the Hilbert space by a factor 2 and 



find no change in the results presented below. In all 
cases the width of the sample L is taken large enough to 
have the cyclotron radius Vc smaller than L/2. The right 
and left edge states are then well separated in real space. 
For yo ~ L/2 the states are equal to the bulk states, 
except for exponential corrections. The energy spectrum 
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FIG. 1: The energy spectrum in units of meV versus the 
guiding center yo without (a) and with (b) the Rashba SOC, 
respectively, for the 2D hole systems. In both figures, the ef- 
fective mass m* —0.27ms, the system size L—600 nm, the hole 
number density Nh=O.Ol nm~^, the magnetic field B—3T, 
and the p-factor (/s=— 0.44. The Rashba SOC strength in 
(b) is set as a=40meVnm^. Clearly from (b) one can see 
that the introduced Rashba SOC mixes the the eigenstates of 
the spin operator Sz- 



the free energy is given by 

F{B, T)^^,M-^^I^' dyo E l^^l + } 

= mA^-^ j dED{E,B)\n{l + eP^^'-''^}, (8) 

where (5=1 /kBT, N^=eB /he, and /i is the chemical po- 
tential. Note that we have defined in the above equation 
the density of states (DOS) per area 

N 

D{E,B) = ^J2 <iyoS{E-E,,^M)- (9) 

n,s -^0 

The explicit inclusion of the DOS in the expression can 
be utilized to take into account the impurity effect, which 
broadens the Landau levels into Gaussian or Lorentzian 
in shape. For simplicity we do not consider the broad- 
ening effect in this paper. In the absence of edge states, 
the Landau Levels En^siuo) a-rc uniform in space and thus 
Eq. (HI) reduces to 



reproduces the bulk results without edge states, which is 
given by the follows: 



^ = ("-l) + ^(l-ff), (" = 0,1,2) (6) 

^ = {n-l) + y^l^n {n - 1) (n - 2) + {3g - 3)^ 

{n > 3) (7) 

with s=±l. The dimensionless parameter is defined by 
7 = 2 (am* /h?) ^^/2eB/hc\2^. As yo approaches the 
sample edge, the effect of the confining potential becomes 
important and it generates the fc-dependent dispersion of 
the energy levels, which has profound effects on magneto- 
transport and magnetization properties. Fig[l] plots the 
the energy spectrum of the Hamiltonian ([T]) as a function 
of the cyclotron center yo- 



III. THE MAGNETIZATION OF THE 2D HOLES 

Now we investigate the magnetization of the 2D holes 
in the III-V semiconductor wells. According to the 
thermodynamics and the statistical physics, one easily 
obtains that the magnetization density is the deriva- 
tive of the Helmholtz free energy density with respect 
to B at fixed electron density J\f and temperature T, 
M = - (dF/dB) l^r^T- For the present 2D hole model. 



T) = - ^ 5] In { 1 + e^e^-^-^) } . (10) 

The B dependent chemical potential fi is connected to 
the experimentally accessible electron density Af via the 
local DOS. In the clean sample limit this is written as 

where f„s (yo)=— jts — ^^ — \ — is the Fermi distribution 

for the spin-split Landau levels i?„,s(yo)- From Eq. © 
the magnetization density becomes 

n.s \^ ^ 

= A/(")+M(i). (12) 

One can see that the magnetization consists of two parts. 
The first part M*^"^ is the conventional contribution from 
the B dependence of the Landau levels and thus is de- 
noted as a paramagnetic response. The second part il/'-'^' 
comes from the B dependence of the level degeneracy 
factor N,y, thus describing the effect of the variation of 
the density of states upon the magnetic field. Obviously, 
AI^^^ is negative while M*^^) is positive, the net result is 
an oscillation of the total magnetization M between the 
negative and positive values as a function of B. At zero 
temperature, the expression for M reduces to a sum over 
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FIG. 2: The dHvA oscillations of the magnetization (per hole) 
of 2D hole system with different system sizes. The parame- 
ters are chosen as follows: the effective mass m*=0.27me, the 
hole number density Nh=0.01 nm~^, the g-factor 0.44, 
and the temperature T = 0.5K. The Rashba SOC here is 
neglected. 



all occupied Landau levels: 



occ f 
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L dB 



(13) 



where the sum runs over all occupied states and /xq is the 
zero-temperature chemical potential (Fermi energy). 

To clearly see the influence of the edge-states and 
Rashba SOC on the magnetization, let us begin with the 
conventional result for the bulk 2D holes without SOC 
and edge-state effects. In this case, the magnetization 
(per hole) m (see the solid line in Fig. ^ displays the 
well-known sawtooth behavior with varying the magnetic 
field B. One fact revealed in Fig. [2] is that the inclusion 
of the Zecman splitting in the Landau levels will result in 
many weak peaks appearing among the dHvA oscillation 
modes of the physical quantities, the chemical potential 
/i and the magnetization (per hole) m, at very low tem- 
perature (here T = 0.2K). These weak peaks have been 
observed recently by Schaapman et when they mea- 
sured the magnetization of a dual-subband 2D electron 
gas, confined in a GaAs/AlGaAs heterojection, and by 
Zhu et al^ when they measured the magnetization of 
high-mobility 2D electron gas. From Fig. [2] one obtain 
that these weak peaks can also be observed in 2D hole 
system. These weak peaks will disappear when the tem- 
perature turns sufficiently high. 

Then we investigate the edge-state effects on the mag- 
netization. Fig. [2] (the dashed, dotted, and dash-dotted 
lines) also shows the influence of the edge states on the os- 
cillations of chemical potential and magnetization (dHvA 
oscillations) with magnetic field. From Fig. [5] one can 



easily observe a prominent feature brought by the edge 
states, which is that the center of the dHvA oscillations 
is now dependent on the magnetic field. In particular, 
for the field less than 1 Tesla, the oscillatory magnetiza- 
tion is always positive in sign. Another feature shown in 
Fig. [2] is that the oscillation amplitude decreases with de- 
creasing the sample size. As is known, the origin of dHvA 
oscillations is the degeneracy of Landau levels. The edge 
states with dispersion lead to edge current, which not 
only is crucial for the quantum Hall effects, but also very 
important for the magnetization The dispersion of 
the edge states partially lift the degeneracy of the Landau 
levels. Thus the edge states tend to destroy the dHvA 
oscillations. Therefore it leads to the decreasing of os- 
cillation amplitude as shown in Fig. [21 The upshift of 
the center of dHvA oscillations may be understood as 
follows: For the effects from the edge states, what really 
matters is the ratio of two important length scales: the 
magnetic length li, and the system size L. The decreasing 
of L is equivalent to the increasing of lb, i.e., decreasing 
of B or (Jc- From Eq. p^ . one can see that with de- 
creasing of L, the second term overcomes the first term 
and leads to the upshift of the center of dHvA oscilla- 
tions. The smaller the system size is, the more profound 
effects the edge states lead to, as shown in Fig. [5]for both 
the center and amplitude of the dHvA oscillations. Fig. 
|3{b) shows quantitatively the system size dependence of 
the shift of the oscillation center. It has the dependence 
1/L. Roughly, the contribution of the edge states is pro- 
portional to the number of edge states (as also seen from 
Eg. p^ ). which is proportional to vTc/L, where the cy- 
clotron radius rc = y/vh [l^, with the number of the 
occupied Landau levels v ~ Thus the center of 

dHvA oscillations is proportional to l^/L = 1/B^L. The 
B and L dependence is clearly seen in Figs. EJa) and (b). 
To see more explicitly the contribution from edge states 
and bulk states, we plot the total magnetization and the 
contribution from bulk states in Fig. HJa). The contri- 
bution from the edge states is obtained from Eq. by 
summing over terms from edge states, with \yo\ < rc or 
\L — yo\ < Tc- The rest contribution is from bulk states. 
There is no upshift of the magnetization oscillation cen- 
ter for the part from bulk states. It shows explicitly that 
the upshift of the center of dHvA oscillations is due to 
the existence of edge states. Fig. IH^b) shows the depen- 
dence of edge states contribution on the magnetic field. 
The contribution from edge states increases as decreasing 
the magnetic field, or equivalently decreasing the sample 
size as one expects. 

When the Rashba SOC is introduced, there is a en- 
ergy competition between the Zeeman coupling and the 
Rashba SOC. Also, due to the entanglement between the 
orbital and spin degrees of freedom, it is difficult to dis- 
tinguish their separate contributions to the total magne- 
tization. These factors make the physical picture of the 
dHvA oscillations to change fundamentally, as shown in 
Fig. [5] for magnetization (per hole) m as functions of B. 
One can see from these two figures that the Rashba SOC 
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FIG. 3: (a) The B dependence of the center of dHvA oscil- 
lations. The system size i=600nm. (b) The dependence of 
the center of dHvA oscillation on the size of the sample. The 
magnetic field B is chosen around 2.0T. In both figures the 
Rashba SOC is neglected. The other parameters are same as 
those in Fig. (2] 




FIG. 4: (a) Gray curve: The bulk contribution of the magne- 
tization m per hole (in units of /is); Black curve: The total 
magnetization per hole m (in units of ^s). (b) The edge con- 
tribution of the magnetization m per hole (in units of fis)- 
In both figures, the Rashba SOC a = and L — 600nm. The 
other parameters are same as those in Fig. (2] 

has no visible influence on the magnetic oscillations of the 
quantity m at large values of B, where the Zceman and 
spin-orbit coupling splitting arc small compared to the 
Landau level splitting. At low magnetic field, however, 
the Rashba SOC modulation of the magnetic oscillations 
becomes obvious, which can be clearly seen by the en- 
larged plots of n and m in the inset in Fig. \5\ for B 
between 1.2T and 1.4T. For comparison, we also re-plot 
in Fig. [S] the cases without Rashba SOC. One can sec 
from these two figures that the SOC brings about two 
new features at low magnetic field: (i) The sawtoothlike 
oscillating structure is inversed, i.e., the location of peaks 
in fi and m with SOC correspond to the valleys without 
SOC. This inversion is due to the different Landau levels 
in the two cases, (ii) The oscillation mode is prominently 
modulated by SOC and a beating pattern appears. This 
beating behavior in the oscillations are due to the fact 
that the Landau levels are now unequally spaced due to 
the presence of SOC. 

Before ending this paper, let us briefly discuss the 
temperature effect. It is well known that the infinite 
temperature will blur the measured physical quantities, 
which including the dHvA oscillations of the magnetiza- 




FIG. 5: The dHvA oscillations of the magnetization (per hole) 
of 2D hole system as a function of external magnetic field B 
with (dashed line, the Rashba SOC strength a — 20 meV nm^) 
and without (solid line) Rashba SOC coupling at T = 0.2K. 
Other parameters are same as those in Fig. [J] The inset is the 
enlarged magnetization oscillations between 1.2T and 1.4T. 
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FIG. 6: The dHvA oscillations of the magnetization (per hole) 
of 2D hole system as a function of the external magnetic field 
B at different temperatures. The Rashba SOC strength is set 
as Q = 20 meV nm'' . The other parameters are same as those 
in Fig. H 

tion. To illustrate the temperature effect, we plot in FiglS] 
the magnetizations of the holes at different temperatures 
T = 0.5, 1.0, 2.0, and 8.0 K, respectively. Obviously with 
the temperature increases, the amplitude decreases and 
the details of the oscillation tend to be blurred. At the 
temperature T = 8. OK, the dHvA oscillations vanish and 
all the details are smeared out. However, the edge-state 
effect, which depends on the sample size and the external 
magnetic field, does not change with the temperature. 

IV. SUMMARY 

In conclusion, we have theoretically investigated the 
dHvA oscillations of the magnetization in the 2D heavy 
hole systems with boundaries in a perpendicular mag- 
netic filed. Especially, we focus on the edge-state effect 
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and the influence of the Rashba SOC on the magnetiza- 
tion. The results show that the effect of the edge states 
and Rashba SOC on the magnetization is less important 
when the magnitude of the magnetic field B is so large. 
However, with B decreasing, the effect becomes evident. 
First, the dHvA oscillation center of the magnetization 
will shift when the edge-state effect is considered. The 
contribution of the edge states to the total magnetiza- 
tion has a roughly linear relation with 1/B'^ and 1/L. 
When the Rashba SOC introduced, the sawtoothlike os- 
cillating structure is inversed and a beating pattern ap- 
pears. The phenomena caused by the Rashba SOC can 
be blurred when the temperature turns high. However, 
the edge-state effect do not change with the temperature 



changing. 
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